Wavelength modulation spectroscopy method and system

ABSTRACT

A method and system for measuring the concentration of a gas component in a measuring gas a provided. The wavelength of a light source is modulated with a modulation signal at a modulation frequency, while the wavelength is swept over an interaction feature of a sample. The intensity of the light source is further modulated at a wavelength outside the interaction feature with a burst signal, where an N-th harmonic of the burst frequency coincides with an M-th harmonic of the modulation frequency. The light is passed to the sample and thereafter to a detector. The detector output is demodulated at the M-th harmonic, and the demodulated detector output is normalized by calculating the ratio between a demodulated detector output portion derived from the light modulated with the modulation signal and another demodulated detector output portion derived from the light modulated with the burst signal.

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims priority of European Patent Office application No. 0702493.7 EP filed Dec. 21, 2007, which is incorporated by reference herein in its entirety.

FIELD OF INVENTION

The invention relates to a wavelength modulation spectroscopy method. It further relates to a wavelength modulation spectroscopy system.

BACKGROUND OF INVENTION

In wavelength modulation spectroscopy (WMS) the wavelength of the light of a tunable light source, such as a diode laser, is modulated with a frequency f₀, while the wavelength is swept over a molecular absorption line of a gas component of interest in a gas sample. As the light propagates along a measurement path through the gas sample, wavelength dependent absorption converts some of the wavelength modulation into an amplitude modulation of the light. Thus, the light will have an overtone spectrum generated by the absorption, the harmonic content of the spectrum being dependent on the width and shape of the molecular absorption line in the gas and the etalons in the optical path of the measuring system. When the light then impinges onto a measuring detector, for example a photodiode, the detector output contains AC components at the modulation frequency f₀ and its higher harmonics Mf₀ (M=2, 3, 4, etc.). Demodulating the detector output at one of said higher harmonics, preferably at 2f₀, shifts the measurement from frequencies near DC, where the light source is noisy, into a higher frequency range, where the noise is lower, thus improving the measurement sensitivity.

SUMMARY OF INVENTION

In order to measure absolute gas concentrations, a suitable normalization method is needed for compensating for general fluctuations in the emitted light intensity and non-gas related transmission in the optical path of the measuring system. For example, in in-situ measurements of trace gases in combustion environments where varying dust loads, high temperature, gas turbulences etc. modulate the light in the kHz range, it is important that the normalization is not distorted by the rapidly changing transmission and turbulences in the measurement path.

Light which propagates through weakly absorbing gases is attenuated exponentially according to the Beer-Lambert law:

$\begin{matrix} {{{I(v)} = {{I_{L}(v)}{T \cdot {\exp\left\lbrack {- {\sum\limits_{i}{c_{i}{\alpha_{i}(v)}L}}} \right\rbrack}}}},} & \left( {{Equation}\mspace{14mu} 1} \right) \end{matrix}$

where I is the intensity of the light after passing through the measurement path, I_(L) is the intensity of the light emitted from the light source, T is a transmission factor over the measurement path, which transmission factor stands for the wavelength independent transmission including optical losses, α_(i) is the absorption coefficient of a gaseous species i with the concentration c_(i), and L is the length of the measurement path. The absorption coefficient α_(i) is dependent on the light frequency v (or the wavelength). For small optical absorption, Equation 1 reduces to:

$\begin{matrix} {{I(v)} = {{I_{L}(v)}{T \cdot {{\exp\left\lbrack {1 - {\sum\limits_{i}{c_{i}{\alpha_{i}(v)}L}}} \right\rbrack}.}}}} & \left( {{Equation}\mspace{14mu} 2} \right) \end{matrix}$

As mentioned above, wavelength modulation spectroscopy utilizes a rapid modulation of the emitted light with a frequency f₀, while the wavelength is swept over a molecular absorption line of a gas component of interest in the gas sample. The light impinging on the detector can then be written as:

I(ν)=I _(L)(ν)T(1−α₀χ(ν)c ₀ L)=I _(BG)(ν)+I _(AS)(ν)  (Equation 3)

[Applied Optics, Vol. 38, No. 27, pp. 5803-5815 (September 1999)], where α₀ and χ represent the intensity and the peak-normalized shape of the molecular absorption line of interest. I_(BG) and I_(AS) are the background and analytical light portions, respectively, and are defined as:

I _(BG)(ν)=I _(L)(ν)T  (Equation 4)

I _(AS)(ν)=−I _(L)(ν)Tα ₀χ(ν)c ₀ L  (Equation 5).

As can be seen from Equations 4 and 5, independent measuring of the non-gas related transmission I_(L)T is needed to measure absolute gas concentrations.

The most straight forward method to measure the non-gas related transmission I_(L)T is to use a direct detection. The wavelength of the light is swept by a triangular or sawtooth waveform over the absorption line of the gas component to be measured wherein the beginning and the end of the scan are well separated from the absorption peak. The measuring detector output is compared with the signal from a monitor detector which directly monitors the output intensity of the light source. The direct detection channel then detects the large triangular scan as a measure of the transmitted optical power. The scan also includes a period where the light source is turned off in order to provide an accurate zero irradiance reference. [Applied Optics, Vol. 38, No. 36, pp. 7342-7354 (December 1999) and Applied Optics, Vol. 44, No. 1, pp. 91-102 (January 2005)].

In wavelength modulation spectroscopy a combination of wavelength modulation and direct detection can be used [Applied Optics, Vol. 38, No. 21, pp. 4609-4622 (July 1999)]. This technique is mostly developed for atmospheric monitoring; to be used in harsh industrial environment, the modulation rate has to be increased in order to place the signal energy above that of the turbulent measuring medium.

In wavelength modulation spectroscopy an indirect measure of the non-gas related optical transmission can be obtained by the use of the wavelength modulation signal f₀, which makes it necessary to introduce a separate detection channel for the fundamental frequency [U.S. Pat. No. 5,173,749]. An intentionally injected pilot tone at a higher harmonic Mf₀ of said wavelength modulation signal [U.S. Pat. No. 7,116,422] avoids the use of such a separate electronic channel. A drawback of this method, however, is that the received pilot tone amplitude gives only information about the transmission factor T rather than the detected non-gas related light intensity I_(L)T, thus I_(L) has to be measured separately, e.g. by division with a reference cell signal [U.S. Pat. No. 5,173,749], which introduces the necessity of an additional optical channel. Therefore, in order to obtain I_(L)T directly, the modulation of the light source should also include turning off the emitted light entirely.

Therefore, the invention seeks to provide a wavelength modulation spectroscopy method and system, which effectively compensate variations in the emitted light intensity and in the non-gas related transmission of the measurement path.

According to the invention this is achieved by the method and the system defined in the independent claims.

Preferred embodiments of the method and the system according to the invention are specified in the remaining claims.

According to the present invention normalization is based on a burst signal, the frequency of which lies above that of the turbulences and flame spectra in the measurement path. The burst signal waveform is optimized to maximize the intensity modulation effect while its amplitude is chosen to allow periodic interruption of the laser emission. Moreover, the burst frequency is chosen so that a suitable overtone can be detected by the same signal chain or channel as that of the analytical signal portion thereby using preferably a down sampling scheme.

BRIEF DESCRIPTION OF THE DRAWINGS

The present invention will be now described by way of example and with reference to the accompanying drawings, in which:

FIG. 1 shows an example for the light modulation with a burst and a sinusoidal signal;

FIG. 2 shows the Fourier spectra of the detected burst and analytical signal portions;

FIG. 3 illustrates the down sampling procedure; and

FIG. 4 is a schematic block diagram of the system in accordance with the invention.

DETAILED DESCRIPTION OF INVENTION

FIG. 1 shows an example of the modulation of the intensity I_(L) of the light emitted from a light source, preferably a diode laser. The modulation periodically alternates between a burst 1, 1′ with a burst frequency f_(n) and a triangular or sawtooth sweep function 2 with a superimposed sinusoidal modulation 3 at a modulation frequency f₀. The purpose of the sweep function 2 is to allow light wavelength scan across an absorption line of a gas component of interest. It is advantageous, although not necessary, to apply the burst 1, 1′ with different amplitudes before and after the sweep function 2 to allow measure of the eventual changes in the optical power vs. injection current characteristics of the diode laser.

The waveform and the frequency f_(n) of the burst 1, 1′ are preferably chosen to allow a settling time of the diode laser. The optimum choice is a square wave. Furthermore, the burst frequency f_(n) is arranged so that an N-th harmonic Nf_(n) of the waveform coincides with an M-th harmonic Mf₀ of the modulation frequency f₀, at which M-th harmonic Mf₀ the analytical light portion I_(AS) (cf. Equation 5) is detected. Thus, for the purpose of normalization, a suitable overtone can be detected by the same signal chain or channel as that for evaluation the analytical signal portion. The choice of harmonics M and N is also made to allow for an optimum relation in amplitude between the normalization signal and the analytical signal, thus allowing an optimum dynamic range of the single signal chain.

The M-th harmonic of the detected analytical light portion (cf. Equation 5) can be written as:

I _(AS,Mf) _(o) (t)=α₀χ_(M) I _(L) Tc _(i) L cos(2πMf ₀ t)  (Equation 6).

The detected light intensity as a result of the square wave burst modulation shown in FIG. 1 can be then written as:

$\begin{matrix} {{{I_{burst}(t)} = {I_{L}T\frac{2}{\pi}\begin{pmatrix} {\frac{\pi}{4} +} \\ {\sum\limits_{{k = 1},3,5,\ldots}^{K}{\left( {- 1} \right)^{{({k + 1})}/2}\frac{1}{k}{\cos \left( {2\pi \; {kf}_{n}t} \right)}}} \end{pmatrix}}},} & \left( {{Equation}\mspace{14mu} 7} \right) \end{matrix}$

where K depends on the bandwidths of the laser driver and the detector.

As the burst frequency f_(n) is arranged so that an N-th harmonic Nf_(n) of the burst frequency f_(n) coincides with an M-th harmonic Mf₀ of the modulation frequency f₀, the following relation is given:

$\begin{matrix} {{f_{0} = {\frac{N}{M\;}f_{n}}},} & \left( {{Equation}\mspace{14mu} 8} \right) \end{matrix}$

where N=1, 3, 5, . . . .

If another waveform configuration is used when the modulation and burst signals simultaneously, a further aspect when determining the relation between these two frequencies f₀ and f_(n) is to avoid distortion due to overlap between frequency components of the measurement and burst signal. In this case the bandwidth B of the measured signal at Mf₀ has to fulfill the following relation:

B<2f_(n)  (Equation 9).

By inserting Equation 8 in Equations 6 and 7, respectively, one obtains:

$\begin{matrix} {{{I_{{AS},{Nf}_{n}}(t)} = {\alpha_{0}\chi_{M}I_{L}{Tc}_{i}L\; {\cos \left( {2\pi \; {Nf}_{n}t} \right)}}}{and}} & \left( {{Equation}\mspace{14mu} 10} \right) \\ {{I_{{burst},{Nf}_{n}}(t)} = {I_{L}T\frac{2}{N\; \pi}\left( {- 1} \right)^{{({N + 1})}/2}{{\cos \left( {2\pi \; {Nf}_{n}t} \right)}.}}} & \left( {{Equation}\mspace{14mu} 11} \right) \end{matrix}$

FIG. 2 shows the Fourier spectra of the detected burst and analytical signal portions I_(burst)(f) and I_(AS)(f) for N=7 and M=2. The Nf₀ frequency components can be filtered and amplified before downsampling.

By performing downsampling at a sampling frequency F_(s)=(N+1)f_(n), Equation 10 can be written as:

$\begin{matrix} {{{I_{{AS},{Nf}_{n}}(n)} = {\alpha_{0}\chi_{M}I_{L}{Tc}_{i}{\cos \left( {2\pi \; n\; \frac{N}{N + 1}} \right)}}},} & \left( {{Equation}\mspace{14mu} 12} \right) \end{matrix}$

where n is a sample number. Similarly, Equation 11 becomes:

$\begin{matrix} {{I_{{burst},{Nf}_{n}}(n)} = {I_{L}T\; \frac{2}{N\; \pi}\left( {- 1} \right)^{{({N + 1})}/2}{{\cos \left( {2\pi \; n\; \frac{N}{N + 1}} \right)}.}}} & \left( {{Equation}\mspace{14mu} 13} \right) \end{matrix}$

Since N/(N+1)>½, aliasing takes place. Reconstruction of the discrete signals given by Equations 12 and 13 gives:

$\begin{matrix} {{{\cos \left( {2\pi \; n\; \frac{N}{N + 1}} \right)}\overset{{f_{n}t} = {n/{({N + 1})}}}{}{\cos \left( {2\pi \; f_{n}t} \right)}}.} & \left( {{Equation}\mspace{14mu} 14} \right) \end{matrix}$

Thus, by performing down sampling at a sampling frequency F_(s)=(N+1)f_(n), the M-th harmonic of the detected analytical light portion I_(AS,Nfn) and the N-th harmonic I_(burst,Nfn) of the burst are both converted down to f_(n) due to the aliasing effect. This effect is shown in FIG. 3.

Combination of Equations 12 and 13 yields the following formula for gas concentration:

$\begin{matrix} {c_{i} = {\frac{2\left( {- 1} \right)^{{({N + 1})}/2}}{\alpha_{0}\chi_{M}\; L\; N\; \pi} \cdot {\frac{I_{{AS},f_{n}}}{I_{{burst},f_{n}}}.}}} & \left( {{Equation}\mspace{14mu} 15} \right) \end{matrix}$

As can be seen, the concentration c_(i) is no longer dependent on the non-gas related optical transmission I_(L)T.

The above method is especially advantageous to utilize an audio analog-to-digital converter with a sampling frequency F_(s)=192 kHz. This avoids the necessity of an extra downsampling stage. The burst frequency is then f_(n)=24 kHz, while the modulation frequency is f₀=84 kHz. The 7f_(n) burst and 2f₀ analytical signal fall both in a 168 kHz frequency band. Sampling at 192 kHz aliases the 168 kHz band back to 24 kHz where they can be easily processed further.

FIG. 4 shows a wavelength modulation spectroscopy system including a frequency tunable light source 11 in form of a diode laser for generating light 12 in form of a laser beam and of intensity I_(L) which is passed along a single optical path through a measuring volume 13 to a detector 14 for generating an output 15 indicative of the received light intensity I. The measuring volume 13, which can be a sample cell or, in case of in-situ process measurements, a gas-leading pipe, furnace, funnel or the like, contains a measuring gas (sample) 16, in which the concentration c_(i) of a specific gas component i is to be measured. The modulation of the diode laser 11 is switched by means of a switch 17 between the sweep signal 2 with the added modulation signal 3 of the frequency f₀, provided by a waveform generator 18, and the burst signal 1, 1′ turning on and off the diode laser 11 at frequency f_(n), provided by a burst generator 19. The frequencies f₀ and f_(n) are related such that Mf₀=Nf_(n), where N is an integer corresponding to a suitable harmonic of the burst signal 1, 1′ and M is the harmonic of the modulation frequency f₀ where detection of the absorption in the measuring volume 13 will be made. The generated laser light 12 is passed through the measuring volume 13 and picked up by the detector diode 14. The output 15 of the detector 14 is filtered through a band-pass filter 20 with a centre frequency Mf₀ and then converted to digital format in an analog-to-digital converter 21 running at a sampling frequency F_(s)=(N+1)f_(n) hence causing both the M-th harmonic of f₀ and the N-th harmonic of f_(n) to be folded or aliased back to frequency f_(n). The down-converted detector output is then processed by the digital signal processing unit 22 for calculating the concentration c_(i) of the specific gas component i to be measured. Due to a synchronisation unit 21 the signal processing unit 22 can separate the parts of the detector output related to modulation generated by the burst generator 19 at the burst frequency f_(n) from those parts related to modulation by the waveform generator 18 at the modulation frequency f₀. 

1.-12. (canceled)
 13. A wavelength modulation spectroscopy method comprising: periodically sweeping the wavelength of a light source over an interaction feature of a sample according to a sweep function; modulating the wavelength of the light source with a modulation signal at a modulation frequency, while the wavelength is swept over the interaction feature; further periodically modulating the intensity of the light source at a wavelength outside the interaction feature with a burst signal, where an N-th harmonic of the burst frequency coincides with an M-th harmonic of the modulation frequency; passing the light of the light source to the sample for interacting and thereafter to a detector; demodulating the detector output at the M-th harmonic of the modulation frequency; and normalizing the demodulated detector output by calculating the ratio between a demodulated detector output portion derived from the light modulated with the modulation signal and another demodulated detector output portion derived from the light modulated with the burst signal.
 14. The method of claim 13, wherein before the step of normalizing, the demodulated detector output is down sampled at a sampling frequency equal to the (N+1)-fold of the burst frequency.
 15. The method of claim 13, wherein N=7 and M=2.
 16. The method of claim 13, wherein the modulation signal comprises a sinusoidal.
 17. The method of claim 16, wherein the burst signal comprises a square wave.
 18. The method of claim 13, wherein the burst signal comprises a square wave.
 19. The method of claim 16, wherein the sweep function comprises a sawtooth, and wherein the amplitude of the burst signal is different before and after the sweep.
 20. A wavelength modulation spectroscopy system, comprising: a wavelength tunable light source; a first modulator that periodically sweeps the wavelength of the light source over an interaction feature of a sample according to a sweep function and modulates the wavelength of the light source with a modulation signal at a modulation frequency, while the wavelength is swept over the interaction feature; a second modulator that periodically modulates the intensity of the light source at a wavelength outside the interaction feature with a burst signal, where an N-th harmonic of the burst frequency coincides with an M-th harmonic of the modulation frequency; a detector that detects the light of the light source after interaction with a sample and producing a detector output; a demodulator that processes the detector output at the M-th harmonic of the modulation frequency, and an evaluator that normalizes the demodulated detector output configured to calculate the ratio between a demodulated detector output portion derived from the light modulated with the modulation signal and another demodulated detector output portion derived from the light modulated with the burst signal.
 21. The system of claims 20, further comprising between the detector and the evaluator an analog-to-digital converter configured to downsample the demodulated detector output at a sampling frequency equal to the (N+1)-fold of the burst frequency.
 22. The system of claims 20, wherein N=7 and M=2.
 23. The system of claims 20, wherein the modulation signal comprises a sinusoidal.
 24. The system of claims 20, wherein the burst signal comprises a square wave.
 25. The system of claims 20, wherein the sweep function comprises a sawtooth, and wherein the amplitude of the burst signal is different before and after the sweep.
 26. The system of claims 21, wherein N=7 and M=2.
 27. The system of claims 26, wherein the modulation signal comprises a sinusoidal.
 28. The system of claims 27, wherein the burst signal comprises a square wave.
 29. The system of claims 28, wherein the sweep function comprises a sawtooth, and wherein the amplitude of the burst signal is different before and after the sweep. 